require(rbacon)
#require(devtools)
#load_all("~/Dropbox/devsoftware/rbacon")
fry <- tempdir()
Bacon(run=FALSE, ask=F, coredir=fry, ssize=500, BCAD=F)

agedepth()
## Calculating age ranges...
## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm
## 100% of the dates overlap with the age-depth model (95% ranges)
## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92
agedepth(rotate.axes=TRUE)
## Calculating age ranges...
## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm
## 100% of the dates overlap with the age-depth model (95% ranges)
## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92
agedepth(BCAD=TRUE)
## Calculating age ranges...
## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm
## 100% of the dates overlap with the age-depth model (95% ranges)
## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92
agedepth(rotate.axes=TRUE, BCAD=TRUE)
## Calculating age ranges...
## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm
## 100% of the dates overlap with the age-depth model (95% ranges)
## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92
agedepth(rev.d=TRUE)
## Calculating age ranges...
## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm
## 100% of the dates overlap with the age-depth model (95% ranges)
## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92
agedepth(rev.age=TRUE)
## Calculating age ranges...
## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm
## 100% of the dates overlap with the age-depth model (95% ranges)
## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92
agedepth(rev.d=TRUE, rev.age=TRUE)
## Calculating age ranges...
## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm
## 100% of the dates overlap with the age-depth model (95% ranges)
## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92
agedepth(rotate.axes=TRUE, rev.d=TRUE, rev.age=TRUE)
## Calculating age ranges...
## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm
## 100% of the dates overlap with the age-depth model (95% ranges)
## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92
agedepth(BCAD=T)
## Calculating age ranges...
## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm
## 100% of the dates overlap with the age-depth model (95% ranges)
## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92
agedepth(BCAD=T, rev.age=T)
## Calculating age ranges...
## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm
## 100% of the dates overlap with the age-depth model (95% ranges)
## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92
agedepth(BCAD=F)
## Calculating age ranges...
## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 206 yr, min. 115 yr at 73.5 cm, max. 341 yr at 25.5 cm
## 100% of the dates overlap with the age-depth model (95% ranges)
## Posteriors: accrate mean 22.27, shape 1.74, memory mean 0.53, strength 9.92
proxy.ghost()
## Calculating histograms
## 

layout(1)
accrate.depth.ghost()

accrate.depth.ghost(rotate.axes=T)

accrate.depth.ghost(rev.d=T)

accrate.depth.ghost(rev.acc=T)

accrate.depth.ghost(rev.d=T, rev.acc=T)

accrate.depth.ghost(rev.d=T, rev.acc=T, rotate.axes=T)

accrate.age.ghost()
## 

accrate.age.ghost(rotate.axes=T)
## 

accrate.age.ghost(rev.age=T)
## 

accrate.age.ghost(rev.acc=T)
## 

accrate.age.ghost(BCAD=T)
## 

accrate.age.ghost(BCAD=T, rev.acc=T, rev.age=T)
## 

accrate.age.ghost(rev.age=T, rev.acc=T)
## 

accrate.age.ghost(rev.age=T, rev.acc=T, rotate.axes=T)
## 

myflux <- cbind(0:100, rnorm(101, 100, 5))
flux.age.ghost(flux=myflux)
## 
## 

flux.age.ghost(rotate.axes=T, flux=myflux)
## 
## 

flux.age.ghost(rev.flux=T, flux=myflux)
## 
## 

flux.age.ghost(rev.age=T, flux=myflux)
## 
## 

Bacon(slump=c(20,21), ssize=500, coredir=fry, ask=F)
##  Will run 2,200,000 iterations and store 500
## Calculating age ranges...
## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 198 yr, min. 98 yr at 74.5 cm, max. 298 yr at 27.5 cm
## 100% of the dates overlap with the age-depth model (95% ranges)
## Posteriors: accrate mean 23.26, shape 1.97, memory mean 0.52, strength 9.25
Bacon("RLGH3", coredir=fry, accept.suggestions=T, ask=F)
##  Using a mix of cal BP and calibrated C-14 dates
## Adapting acc.mean to 50 yr/cm
##  Will run 9,187,500 iterations and store 4,000
## Calculating age ranges...
## 
## Preparing ghost graph...

## 
## Mean 95% confidence ranges 578 yr, min. 57 yr at 0 cm, max. 1835 yr at 220 cm
## Warning! Only 76% of the dates overlap with the age-depth model (95% ranges)
## Posteriors: accrate mean 46.08, shape 0.81, memory mean 0.47, strength 12.25
accrate.age.ghost(BCAD=T)
##